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Surface- induced  optical  bistability  in  coupled  exciton- phonon  systems 
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Abstract 

Nonlinear  optical  responses  of  exciton-phonon-coupling  systems  near  a 
reflecting  surface  are  investigated.  For  a  polydiacetylene  chain  lying 
parallel  to  the  surface,  optical  bistability  induced  by  the  surface  is 
discovered.  It  is  also  found  that  the  threshold  and  contrast  of  the 
bistability  can  be  controlled  by  adjusting  the  distance  of  the  system  from  the 
surface,  and  that  vacuum  fluctuations  are  reduced  whenever  the  bistability 
occurs.  The  nature  of  the  bistability  and  its  origin  are  discussed. 


1990  PACS  Numbers:  42.65. -k,  42. 65. pc,  36.20. -r,  78. 65. He 


Ever  since  the  mid- 1940 ' s ,  it  has  been  known  that  the  spontaneous 

emission  probability  of  a  spin  system  can  be  dramatically  influenced  by  its 

environment.  In  recent  years,  inhibited  spontaneous  emission  has  been 

2  3 

investigated  extensively  in  atomic  systems  and  in  solids.  It  is  found  that 

the  spontaneous  radiation  rate  is  altered  markedly  by  merely  changing  the 

surrounding  material  with  a  different  index  of  refraction.  On  the  other  hand, 

a  great  deal  of  work  has  been  devoted  to  the  study  of  an  adatot:  near  a  tclid 
4  5 

surface.  ’  Optical  properties  are  radically  changed  because  of  the  presence 
of  the  surface,  and  the  change  depends  on  the  structure  of  the  substate 
material.  Such  problems  are  of  fundamental  importance  as  well  as  of  practical 
interest.  Fundamentally,  one  has  to  develop  a  quantum  mechanical  theory  for 
dissipative  systems  to  deal  with  such  nonlinear  interactions  between  the  atom 
and  the  substrate  which  does  not  possess  bulk  symmetry.  From  a  practical 
point  of  view,  such  studies  car  lead  to  applications  in  surface  photochemistry 
and  surface  optics. 

Instead  of  atoms  or  molecules  adsorbed  at  a  solid  surface,  we  consider 
an  exci ton- phonon  coupling  system  such  as  a  one -dimensional  semiconductor  of 
polymers.  In  particular,  we  consider  a  sample  of  polydiacetylene  (PDA) 
located  near  a  metal  surface.  The  system  closely  resembles  the  experimental 
situation^  in  which  PDA  films  are  sandwiched  between  other  materials.  PDA  is 
a  good  candidate  for  such  considerations  because  of  its  large  third-order 

/  3  \  7  g 

susceptibility  x  .  ^nd  because  of  its  small  transmission  loss  a.  Thus  it 

(3) 

possesses  a  fairly  large  ratio  x  /<» .  which  is  usually  a  measure  of  the 

usefulness  of  a  material  employed  in  switching  devices.  Moreover,  PDA  is  more 
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flexible  to  use  in  the  construction  of  waveguides. 

For  a  linear  chain  of  PDA- toluene-sulfonate  (PTS)  placed  at  a  distance  d 
from  and  oriented  parallel  to  the  metal  surface,  we  shall  study  its  response 


to  a  laser  beam  shining  on  it.  Novel  phenomena  have  been  found  recently  in 

the  study  of  the  nonlinear  optical  response  of  PTS  in  a  thin  film  in  the 

transient  regime,  and  in  an  optical  cavity  in  the  steady-state  regime. 

In  this  Letter,  our  purpose  is  to  investigate  the  possibility  of  optical 

multistability  induced  by  the  metal  surface,  even  though  there  is  no  cavity  to 

feedback  the  incident  signal.  The  light- induced  excitons  in  PTS  act  as 

emitting  dipoles.  The  emitted  light  is  reflected  by  the  metal  surface  and 

interacts  with  the  excitons  themselves,  and  changes  their  dynamical  behavior. 

12  lb 

In  addition,  the  excitons  are  also  coupled  to  phonon  modes  in  PTS.  '  Such 


coupling  plays  a  very  important  role  in  the  nonlinear  optical  response  of 


polymers 
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This  complicated  interacting  system  is  treated  by  modeling  the  excitons 

and  phonon  modes  by  damped  oscillators  and  by  assuming  dissipative 
14 

interactions  between  the  exciton  and  the  reflected  field.  Thus  the  problem 
is  described  by  the  nonhermitian  Hamiltonian 


H-H  +  H  +  H  +H^  +  H„  , 
X  p  xp  xf  xR 


where  the  free  exciton  energy  operator  is 


the  free  phonon  energy  is 


“p  -  K )  ''Ih  • 


(3) _ 

on/_ _ 

Availability  Codes 
lAvall  and/or 
01st  I  Speoial 


□  □ 


the  exci ton -phonon  interaction  energy  is 


H  -  >(  )  A.  a^a(b.+b!)  ,  (4) 

xp  '  /  1  1  1 

i 

the  interaction  energy  between  the  exciton  and  driving  field  in  the  rotating- 
wave  approximation  (RWA)  is 

*  *  t 

-  -(fiaLE^e  +  n  ^i'EqS  )  ,  (5) 

and  the  interaction  energy  between  the  exciton  and  surface-reflected  field  in 
the  RWA  is 

"xR  - 

The  notation  is  as  follows.  The  operators  a^(a)  and  b!(b.)  create 

(annihilate)  an  exciton  with  energy  and  a  phonon  in  the  i  -  th  mode  with 

energy  yiu. ,  respectively.  The  corresponding  decay  rates  are  denoted  by  y  and 

7^,  respectively,  and  A^^  is  a  coupling  constant.  The  dipole  moment  is  given 

as  p  -  ^a  with  p  as  its  matrix  element,  and  the  driving  field  is  E  = 

EQexp(iwQt)  with  amplitude  Eq  and  frequency  Wq.  stands  for  the  reflected 

field  at  the  position  of  the  dipole.  It  is  proportional  to  p  and  hence  can  be 

written  as  E„  -  E  a,  where  E  is  a  c -number.  As  we  shall  see  later,  E„  is  in 

general  a  complex  quantity,  so  that  the  interaction  (6)  is  dissipative.  We 

also  note  that  the  momentum  dependence  of  the  exciton  has  been  neglected  in 
10-12 


this  Hamiltonian. 


Following  Dekker's  quantization  procedure  for  dissipative  systems, the 
equation  of  motion  for  the  density  operator  of  the  PTS  absorbate  can  be 
written  as 

P  -  -  ^  [a^[a,H]p]  +  ^  [p[H^a^],al 

i  i 

Since  we  are  interested  in  the  scattered  field  intensity  I  as  a  function  of 

sc 

the  input  optical  field  intensity  the  effects  of  quantum  fluctuations  are 

ne g  1 i gib  1 e  ,  ^ ^  and  a  s emiclass ical  approach  is  applicable.  Thus  it  is 

sufficient  to  use  the  mean  values  a  =  <a>,  n  =  <a^a>  and  B.  =  <b.>.  The 

11 

equation  of  motion  for  these  variables  in  the  Schrodinger  picture  can  be 
obtained  directly  from  Eq.  (7).  In  the  rotating  frame,  we  have 


•  • 


a  -  <a>  =“  Tr(pa) 


=  -  1 


1  Y  * 

(A+w  )  +  ':77  +  7  ]«  +  iO  -  i)  A.(^.+^.)a 

'  s  2  X  s'  /  1 


(8) 


n  -  -(7  +27  )n  -  in  a  +  iOa 

X  s 


(9) 


-  -(ivPi^^i  -  - 


(10) 


where  we  have  defined  the  detuning  A  -  “x'^0  Rabi  frequency  Ci  -  u  Eq- 

The  surface  effect  on  the  exciton  is  reflected  in  the  surf  ace  -  induced 


frequency  shift 
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-  |M|^Re(E^/M)  (11a) 

and  Che  decay  rate 

7^  -  \n\hm(E^/^i)  .  (11b) 


The  reflected  electric  field  in  the  dipole  direction  can  be  found  by 
a  classical  approach.  We  take  the  oscillating  dipole  moment  to  be  located  at 
a  distance  d  from  the  surface  of  a  semi- inf inite  metal.  From  electromagnetic 
field  theory,  we  find  in  a  straightforward  manner  that 


2  ^ 


f 


du  u[Rj_  +  (l-u^)R,l  j/yi-u^  (12) 


when  the  dipole  p  is  parallel  to  the  surface,  and 


1 

E 

r 


du  u^R||/yi-u^ 


(13) 


when  the  dipole  is  perpendicular  to  the  surface.  Here  we  have  defined  = 
w/c  and 


Rii  - 


h^2  '  ^2^1 
^1^2  ^2^1  ’ 


(14a) 


R 


1 


~  ^2 
^1  ^^2 


(14b) 


^2  - 


2  -  1,2 


(14c) 


with  £  as  the  dielectric  constant  and  the  subscripts  1  and  2  labelling  the 
medium  containing  the  PTS  chain  and  the  metal  substrate,  respectively.  For 
definiteness,  we  consider  the  case  where  p  is  parallel  to  the  surface.  It 
then  follows  directly  from  Eqs.  (11)  and  (12)  that 


3  r  2  ^^i"" 

-  7q  Im  J  du  u[(l-u  )R||  +  R^je 


(15a) 


3  r®  2  ^^1^ 

-  -  -  7q  Re  J  du  u[(l-u  )R||  +  R^je 


(15b) 


where 


u> 


Xv3 


^  “  3  ^  ^c  ^ 


2  \ 


(16) 


is  the  decay  rate  in  the  absence  of  the  surface, and  x  -  2dy7^w^/c. 

In  the  steady  state,  the  solution  for  a  can  easily  be  found  from  Eqs 
(8)-(10)  as 


a  -  n/[A+a)^-A  n-i(7g+p^): 


(17) 


where  A  is  defined  by 
P 


i  i 


because  »  7^ . 


(18) 


The  mean  number  of  excitons  is  determined  by  the  equation 


(19) 


-  2(A+w  )A  +  [  (A+tj  )^  +  (7  +■“7  )^]n  -  I.  -  0 
p  s  p  s  's  2  'x^  ^  in 

2 

where  the  input  laser  intensity  is  1^^^  -  |n|  In  the  scattering  region  where 

the  incident  field  vanishes,  it  is  known  that  the  positive  -  frequency  part  of 

18 

the  field  may  be  written  as 


E^'^\?,t) 


(r,t) 


2 

CJ 

X  r 

—  (p  X  r)  X  r  a(t--) 
4»rc  r 


(20) 


where  E^(r,t)  describes  the  free  propagation.  The  mean  value  of  the  scattered 
intensity,  in  the  steady  state  is  thus  given  by 


-  I,  /{[A+w  -A  n(I,  )]^  +  (7  +^7  )^) 
in'  ^  s  p  '  in  ‘  's  2  x 


(21) 


where 


2 

Io(r)  “  [ - ^  (p  X  r)  X  r] 

4jrc  r 


(22) 


The  intensity  is  studied  numerically  as  a  function  of  according 

to  Eq.  (21)  for  a  PTS  chain  parallel  to  the  surface  of  an  ideal  metal  for 

which  -*  ’*>.  As  a  realistic  approximation,  we  have  assumed  that  only  four 
13 

phonon  modes  are  strongly  coupled  to  excitons  in  PTS.  The  parameters  for 
these  modes  are:^^  -  5.16,  A^  -  2.0;  (*>2  -  3.68,  A2  ”  1.66;  ”  2.98, 

A^  ”  0.46;  and  -  2.36,  A^  -  0.48.  For  simplicity  we  have  chosen  ^7^  to  be 


rhe  unit  for  all  quantities  of  the  dimension  t  in  our  calculations.  Optical 

bistability  is  found  when  the  input  field  frequency  matches  the  exciton 

frequency.  Some  of  our  results  for  the  resonant  excitation  (A  -  0)  are 

plotted  ill  Fig.  1,  clearly  demonstrating  that  the  bistability  is  induced  bv 

the  surface.  Our  numerical  calculation  shows  a  strong  dependence  of  the 

optical  bistability  on  the  distance  d  of  PTS  from  the  surface,  especially  when 

X  <  2.5.  As  X  decreases  further,  it  is  observed  that  the  contrast  of  the 

bistability  (the  ratio  of  the  intensity  of  the  upper  branch  to  that  of  the 

lower  branch)  increases  markedly,  and  that  the  threshold  (the  smallest  input 

field  intensity  for  multiple  solutions  of  I  )  decreases  at  the  same  time. 

sc 

Furthermore,  we  have  also  found  that  the  optical  bistability  appears  only  when 
7^  becomes  negative.  In  other  words,  reduced  vacuum  fluctuations  alwavs 
accompany  the  optical  bistability. 

2  19  17 

It  has  been  established  both  experimentally  ’  and  theoretically  that 

the  spontaneous  emission  changes  dramatically  when  the  atomic  system  is  near  a 

reflecting  surface  or  between  two  mirrors.  This  implies  that  the  mode 

structure  of  the  vacuum  field  suffers  significant  modification  because  of  the 

finite  geometry.  A  similar  situation  in  solid-state  physics  and  electronics 

has  also  been  reported,  and  the  spontaneous  electron-hole  recombination  rate 

can  be  controlled  by  adjusting  the  mirror  position  or  by  varying  the  index  of 

3 

refraction  of  the  medium.  What  we  have  found  here  is  that  by  adjusting  the 
environmental  conditions  to  reduce  the  radiation  rate,  one  can  not  only 
produce  low-noise  lasers  from  a  semiconductor  quantum  well  but  also  create 
optical  bistability  from  a  polymer  chain  near  a  reflecting  surface. 

The  origin  of  the  optical  bistability  may  be  understood  by  looking  at 
the  driving  frequency  dependence  of  the  scattered  field  intensity.  For  a 


10 


fixed  driving  field  intensity,  its  frequency  can  be  scanned  over  the  near¬ 
resonance  region.  Thus  the  exciton  number  becomes  a  function  of  oJq,  as  can  be 
seen  from  Eq .  (19).  Consequently,  can  be  regarded  as  a  function  of 

also.  Figure  2  shows  the  surface- induced  optical  bistability  from  the  point 
of  view  of  the  driving  field  frequency  instead  of  its  intensity.  It  is  found 
that  there  are  three  possible  absorption  rates  for  a  certain  small  frequency 
range.  The  scattered  intensity  in  the  higher  branch  corresponds  to  the  higher 
absorption  rate,  and  in  the  lower  branch  to  the  lower  absorption  rate.  Since 
the  absorption  part  of  the  optical  susceptibility  is  proportional  to  the 
scattered  field  intensity  for  fixed  input  intensity  and  distance  x,  the 
bistability  described  above  is  mainly  an  absorptive  one.  In  other  words,  the 
surface- induced  absorption  rate  is  responsible  for  this  particular  optical 
bistability.  This  is  easily  understood  from  energy  considerations,  because  a 
higher  output  intensity  of  the  excitonic  dipole  results  from  more  effective 
absorption  of  the  input  energy,  and  vice  versa. 
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Figure  Captions 

1.  Scattered  field  intensity  I  vs  the  driving  field  intensity  I.  for 

sc  °  -'in 

various  distances  x  when  A  -  0:  (a)  X  -  CO,  (b)  X  -  1.5,  (c)  x  -  1.1, 

(d)  X  -  0.9  and  (e)  x  -  0.85. 

2.  I  as  a  function  of  the  driving  field  frequency  for  a  fixed  I. 

sc  ^  ^  0  in 

0.10:  (a)  X  -  CO  and  (b)  x  -  1.5. 
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